Existence of exotic torus configuration in high-spin excited states of Ca 
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We investigate the possibility of the existence of the exotic torus configuration in the high-spin excited 
states of 40 Ca. We here consider the spin alignments about the symmetry axis. To this end, we use a three- 
dimensional cranked Skyrme Hartree-Fock method and search for stable single-particle configurations. 
We find one stable state with the torus configuration at the total angular momentum J = 60 % and an 
excitation energy of about 170 MeV in all calculations using various Skyrme interactions. The total 
angular momentum J = 60 h consists of aligned 12 nucleons with the orbital angular momenta A = 
+4, +5, and +6 for spin up-down neutrons and protons. The obtained results strongly suggest that a 
macroscopic amount of circulating current breaking the time-reversal symmetry emerges in the high-spin 
excited state of 40 Ca. 

PACS numbers: 21.60 Jz, 21.60.Ev, 27.40.+Z 



The investigations of nuclei rotating extremely fast about 
the symmetry axis provide a good opportunity to severely 
test the fundamental theory of a quantum mechanics. In 
a classical picture for such rotation, the nuclear density at 
the equatorial plane increases due to the strong centrifugal 
force, and the oblate deformations develop with increas- 
ing rotational frequencies [1]. However, such a collective 
rotation about the symmetry axis is quantum-mechanically 
forbidden. Instead, the spins and the orbital angular mo- 
menta of single particles are aligned with the symmetry 
axis building extremely high angular momenta 12,0]. In 
this case, the total angular momentum of the nucleus is just 
the sum of the symmetry-axis components of those of the 
single particles. Bohr and Mottelson pointed out that in 
such a nucleus, a "macroscopic" amount of circulating cur- 
rent breaking the time -reversal symmetries emerges, which 
is a fascinating new form of the nuclear matter yfl. 

A typical example of such phenomena is the high-K 
oblate isomer states in 152 Dy [5|]. In the high angular- 
momentum states above / = 14 H, the observed excited 
states are irregularly distributed around the average yrast 
line, which is calculated by assuming a macroscopic nu- 
clear shape with the strongly oblate deformations (J3 ~ 
-0.3). The emergence of these high-spin isomers around 
the yrast line can be naturally explained by the single- 
particle alignments [3], because the E2 transitions, which 
are a characteristic quantity related to the rotational col- 
lectivity, are strongly forbidden for those nuclei. In addi- 
tion, the y transitions are strongly suppressed due to the 
large difference of the internal structure between the initial 
and final states. Those high- spin isomers around the yrast 
line are thus called "yrast traps". Many experiments have 
been attempted to produce yrast traps with extremely high 
spins [6]. The observed highest angular momentum in this 



category is / = 49 % in 158 Er 0]. 

In this Letter, we investigate the possibility of the exis- 
tence of the torus configuration as an extreme limit of the 
high-spin oblate isomers. In such a limit, more nucleons 
with higher orbital angular momentum are aligned and its 
density is much denser at the equator part. We here con- 
sider a special case for the existence of the density only 
around its equator part. We below show that such exotic 
state can indeed exist at an extremely high angular momen- 
tum in 40 Ca. 

Many theoretical calculations have been performed to 
search for the high spin states with strongly oblate defor- 
mations in a wide range of nuclei. Those studies have 
been mainly focused on the nuclei with the oblate defor- 
mation at around /? 0.3 heavier than 40 Ca [4]. In con- 
trast, Wong discussed the stability of the torus configura- 
tion using the macroscopic-microscopic model in heavy- 
mass systems 1181- 41 Oil . The calculations using the constraint 
Hartree-Fock (-Bogoliubov) models have been also per- 
formed for heavy-mass systems ifTTlfrHl . The calculations 
using a cluster model have been attempted to obtain the 
ring configurations consisting of alpha particles and extra 
neutrons II 1 311 . However, as shown in the present study, 
an essential mechanism for the stability of such configu- 
ration is the extremely high angular momentum directed to 
the symmetry axis, resulting in the breaking of the time- 
reversal symmetry in the intrinsic state. 

To investigate the possibility of the existence of the torus 
configuration in 40 Ca with high angular momentum, we use 
the three-dimensional Skyrme Hartree-Fock (HF) method 
with a Lagrange multiplier, which is introduced to obtain 
the single-particle states aligned with the symmetry axis. 
To this end, we minimize the HF Hamiltonian, H, with the 
Lagrange multiplier, co. In the present study, we take the 
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FIG. 1. (color online) Total density for (a) the initial condition 
of the HF iterations and (b) the calculated result with tia> = 1.5 
MeV at the 15000 HF iterations. The density is integrated in the 
z direction. The contours correspond to multiple steps of 0.05 
for 2 . The color is normalized by the largest density in each plot. 



z axis as the symmetry axis. We first define the effective 
Hamiltonian, H', given by H' = H - a>J z , where J z de- 
notes the operator for the sum of the z components of the 
total angular-momentum for each single particle, j z , given 

by Jz = HiJz - 1° tn e HF approximation, H' is rewritten 

as H' = Y,Ahi — ^f z ), where hi denotes the Hamiltonian 
for each single particle. The eigenvalue of H' is given by 
(H'j = J]i{( e i ~ X) ~ f ia >Q), where A denotes the Fermi 
energy at a> = and and £2, denote the energy and the z 
component of the total angular momentum in the unit of fi 
for each single particle, respectively. In the present study, 
we search for the stable state using the equivalent cranked 
Skyrme HF equation, 5(h - cdJ z ) = HO], by scan- 
ning a large range of a>. 

Before the HF calculations, we here discuss the shell 
structure of the torus configuration using the radial dis- 
placed harmonic oscillator (RDHO) model |9|]. For the 
torus configuration, not only Q but also the z component of 
the orbital angular momentum, A, are good quantum num- 
bers (Q = A + S, where E denotes the z component of the 
spin values, ±1/2). Two nucleons in each A energetically 
degenerate with the different spin values. At tun = 0, the 
lowest configuration for 40 Ca is A = 0, ±1, ±2, ±3, and ±4 
and the residual two nucleons can occupy any two states 
with A = ±5. At hco + 0, the possible spin aligned config- 
urations are (i) A = 0, ±1, ±2, ±3, ±4 and +5 for the total 
angular momentum J = 20 h[= 5 fix 2 (spin degeneracy) 
X 2 (isospin degeneracy)], (ii) A = 0, ±1, ±2, ±3, +4, +5 
and +6 for J = 60 ft [= 15 % x 2 x 2], and (iii) A = 0, ±1, 
±2, +3, +4, +5, +6 and +7 for J = 100 fi [= 25 fi x 2 x 2]. 

In the self-consistent calculations, the single-particle 
wave functions are described on a Cartesian grid with a grid 
spacing of 1 .0 fm. We take 32 x 32 x 24 grid points for the 
x, y, and z directions, respectively. This was sufficiently ac- 
curate to provide converged configurations. The damped- 
gradient iteration method iflfjl is used, and all derivatives 
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FIG. 2. Convergence behavior of the expectation value of J z in 
the HF calculations versus the number of iterations. The dashed, 
dash-dotted, solid, and dotted lines denote the calculated results 
with %oj = 0.5, 1.0, 1.5, and 2.0 MeV, respectively. 



are calculated using the Fourier transform method. We take 
three different Skyrme forces which all perform well con- 
cerning nuclear bulk properties but differ in details: SLy6 
as a recent fit which includes information on isotopic trends 
and neutron matter II 1711 . and SkI3 and SkI4 as recent fits es- 
pecially for the relativistic isovector structure of the spin- 
orbit force IU8I1 . However, except for the effective mass, 
the bulk parameters (equilibrium energy and density, in- 
compressibility, and symmetry energy) are comparable in 
the all interactions. 

For the initial wave functions, we chose the ring configu- 
ration with 10 a particles placed on the x-y plane, as shown 
in Fig. 1(a). Each a particle is described by the Gaussian 
function with its center placed on z = 0. Using this ini- 
tial condition, we perform the HF iterations with 15000 
times and investigate the convergence of the calculated re- 
sults. Figure 2 shows the convergence behaviors of (j z ^ 
versus the number of the HF iterations with various o/s 
using the SLy6 interaction. We can see that the result cal- 
culated with %oj = 1.5 MeV converges rapidly to J z = 60 %. 
Figure 1(b) shows the density obtained with fia) = 1.5 MeV 
at the 15000th iteration step. The calculated result is in- 
deed the torus configuration. The obtained density distribu- 
tion, p(r, z), can be well fitted by p(r, z) = poe~^ r ~ r ° ) ' +z '^ a ~ , 
where po = 0.13 fm" 3 , ro = 6.07 fm, and cr = 1.61 fm for 
the SLy6 interaction. 

On the other hand, the calculated results with %u> = 0.5 
and 1.0 MeV lead to unstable states. That for ha> = 2.0 
MeV leads to the fission. Although it seems that (-/^)'s 
with fiu = 0.5 and 1.0 MeV converge at the 15000th step, 
those are in fact unstable. In Fig. 2, we can see that these 
states first converge to a quasi-stable state. After that, the 
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FIG. 3. Single-particle energies versus the z component of the 
total angular momentum Q. The solid and open circle denote the 
single-particle states with the positive and the negative parties, 
respectively. The single-particle states with the z component of 
the spin S = +1/2 and -1/2 are connected by the dashed and 
dotted lines, respectively. The solid and the dash-dotted lines 
denote the sloping Fermi energies with Tiu> = 1.2 and 1.7 MeV, 
respectively. The gray area denotes the states below the Fermi 
energy. 



instability of those states increases. In those quasi-stable 
states, J z is 20 fi for fiu = 0.5, and 1.0 MeV and is 100 fi 
for ha> = 2.0 MeV. Later, we will discuss those quasi-stable 
states with J z = 20 and 100 fi. 

We next investigate the region of the Lagrange multipli- 
ers where the torus configurations stabilize, using various 
Skyrme interactions. We find that such region of ha> for the 
SLy6 and SkB interactions extends from 1.2 to 1.7 MeV. 
That for the SkI4 interaction is from 1 .2 to 1 .8 MeV. In such 
regions, we obtain the only one stable torus configuration 
at ./, = 60 fi in all calculations. The obtained excitation en- 
ergies are 170.45, 174.22, and 172.53 MeV for the SLy6, 
SkB, and SkI4 interactions, respectively. The fitted radius 
of the torus configuration is almost the same for all cal- 
culations. The dependence of the calculated result on the 
choice of the interaction is weak. This is mainly because 
that the spin-orbit forces in the radial directions of the inner 
and outer surface of the torus configuration almost cancel 
out each other (see Eq. (4. 14) in Ref. @]). Thus, the RDHO 
mode is a good approximation for the present calculations. 

Figure 3 shows the single-particle spectrum versus Q. 
calculated using the SLy6 interaction. The solid and open 
circles denote the single-particle states with the positive 
and the negative parities, respectively. The single-particle 
states with Z = +1/2 and -1/2 are connected by the dashed 
and the dotted lines, respectively. These energies con- 
nected by each line are well proportional to A 2 , which is 
consistent with that of the RDHO model [see Eq. (4.26) in 
Ref. [9J]]. We here define the "sloping" Fermi energy, A', 





FIG. 4. (color online) Circulating currents of 12 nucleons. The 
outer layer of the torus surface denotes the half of the total den- 
sity. The density plot in the inside of the torus surface denotes 
the calculated density of the aligned 12 single particles. 



given by A' = A + fiioQ. |4|. Then, the occupied states are 
obtained by (//') = - A') with <?, < A'. The solid 

and the dash-dotted lines denote the sloping Fermi ener- 
gies with fiuj = 1.7 and 1 .2 MeV, respectively. The occu- 
pied states are denoted by the gray area. In Fig. 3, we can 
clearly see that the torus configuration can exist at J z = 60 
h in the region from ftu> = 1 .2 to 1.7 MeV. The total angu- 
lar momentum of J z = 60 fi consists of the spin-up-down 
pairs of the aligned single particles with the same A, (£1+ , 
Q_), of (+9/2 + , +7/2 + ) with A = 4, (+11/2", +9/2") with 
A = 5, and (+13/2 + , +ll/2 + ) with A = 6 (see the solid 
ellipses in Fig. 3). 

As discussed in the shell structure of the torus configura- 
tion, the states with J = 20, 60, and 100 % are the pos- 
sible combinations of the spin alignment for the RDHO 
model. For the quasi-stable state with J = 20 h, the cen- 
trifugal force is insufficient to stabilize the torus configura- 
tion against the strong nuclear attractive force. On the other 
hand, to build the state with J = 100 fi, the unbound states 
(e > 0) with D, = +15/2 and +11/2 are occupied, result- 
ing in the instability for the torus configuration. Only for 
J = 60 fi, the torus configuration can be stabilized. Thus, 
the stability of the state with J = 60 fi is rather robust. 
Although we have performed the similar calculations for 
24 Mg and 32 S, we could not find the stable torus states for 
those nuclei. 

As shown in the present study, if the state with the torus 
configuration is formed at J = 60 fi, the macroscopic cir- 
culating current of 12 nucleons with A = +4, +5, and +6 
strongly violates the time -reversal symmetry in the intrin- 
sic state (see Fig. 4). It is interesting to investigate how this 
fascinating new state can be observed in experiments. The 
state with the torus configuration at J = 60 fi would have 
an extremely large magnetic moment (// = 30//^). This 
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would lead to a procession motion under an external mag- 
netic field. 

A question then arises how such a "femto-scale mag- 
net" rotates spontaneously. When the spherical symmetry 
is broken, the collective rotation emerges spontaneously, in 
principle, to restore the broken symmetry. However, it is 
unclear whether such a state built with significant amount 
of circulating current can rotate about the perpendicular di- 
rection to the symmetry axis or not. If not, such a state 
would be an anomalous one, which has not yet been recog- 
nized in the experiments. Even if the state can rotate, the 
rotational band built on such a state would show interesting 
behaviors for their M1/E2 transition strengths and moment 
of inertia, It is intriguing to investigate the extent of 
the difference between the experimental and the classical 
rigid-body moments of inertia (J^]_ = 20.97 ft 2 -MeV" 1 in 
the case of the SLy6 interaction). 

To identify the torus configuration in the experiment, it is 
also important to discuss the competition to the fission de- 
cay channel. The calculation using the macroscopic model 
with the FRLDM2002 parameter set |Q]1 shows that the 
fission barrier for the spheroidal deformations vanishes at 
J = 32 ft and an excitation energy of 49.27 MeV 112X311 . Our 
calculated results for the torus configuration are consider- 
ably higher than those values. However, it would be pos- 
sible that the state with the torus configuration can survive 
against the fission. One important example for such pos- 
sibility is the long-lived K* = 16 + isomer in 178 Hf with a 
half -lives of about 31 y lulll . which is extremely long com- 
pared to other isomers (4 s for the K n = 8" isomer). In 
this connection, it is unclear how the collective path from 
a topologically different torus configuration is connected 
to the region of the spheroidal deformations leading to the 
fission floL fl2h . Wong indeed showed that the torus con- 
figurations are stable from J ~ 57 to 74 % for mass number 

A ~ 50 m. 

In summary, we have suggested the existence of the torus 
configuration in the extremely high- spin excited states 
of 40 Ca using the three-dimensional cranked Skyrme HF 
method. We found only one stable state with the torus 
configuration at J z = 60 fi in all the calculations with any 
Skyrme interactions. The calculated excitation energies of 
this state are 170.45, 174.22, and 172.53 MeV for the SLy6, 
SkI3, and SkI4 interactions, respectively. The obtained re- 
sults are insensitive to the choice of the Skyrme interac- 
tions, because the contribution of the spin-orbit force is 
small in the torus configuration. To build the torus state 
with J = 60 h, the 6 nucleons for D. = +9/2 + and +7/2 + 
with A = 4, £1 = +11/2" and +9/2" with A = 5, and 
£2 = +13/2 + and +ll/2 + with A = 6 are aligned with 
the symmetry axis for both protons and neutrons. We have 
shown that such configuration is built on the major shell 



structure estimated by the RDHO model for the rotating 
torus shape, indicating that the torus state is robustly stable 
and the macroscopic amount of the time -reversal symme- 
try breaking occurs. Although the observation of the torus 
state would be difficult, exploration for such an exotic state 
would provide us valuable information on the new frontier 
of the nuclear matter, which is a big challenge both theo- 
retically and experimentally. 
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